Electronic Signature of the Intermediate Phase in Ge^Sei-^ Glasses: 
A Joint Theoretical and Experimental Study 
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Raman and calorimetric studies on Ge x Sei_ x glasses have provided evidence for the existence of 
the intermediate phase (IP) in chalcogenide and other glasses. Here, we present X-Ray Absorption 
Near Edge Structure (XANES) measurements on germanium selenide glasses in the IP composition 
range, and detect an electronic signature of the IP. A b initio molecular dynamics (MD) based models 
of these glasses are discussed, and an atomistic picture of the IP, based upon the models and available 
experiments is presented. We show that these models reproduce detailed composition-dependent 
structure in the XANES measurements, and otherwise appear to properly represent the structure 
of the GeSe glasses near the IP. 



Few binary materials enable scientific inquiry across 
a continuous range of compositions. This is inconceiv- 
able in crystalline systems, which, even in the most com- 
plex materials, have a finite number of phases available. 
An alternative is to study amorphous materials, where 
glass formation is usually possible for continuous compo- 
sition ranges. An attractive example of such a system is 
Ge^Sei-a;. In sustained study of these and other chalco- 
genides, Boolchand [l[ discovered that various exper- 
imental observables display consistent and anomalous 
behavior through a finite composition window near the 
floppy to rigid transition of Thorpe (also associated 
with the work of Phillips H). This window has now been 
observed in many glasses, and has come to be known as 
the Intermediate Phase (IP). The name originates in the 
idea that within the window, the system is sandwiched 
between a "floppy" and stressed rigid phase. Boolchand 
first observed the IP in Raman studies of tetrahedral 
breathing modes in Ge^Sei-^ glasses for x G (0.20, 0.25). 
129 I Mossbauer measurements reveal [5| [6[ that the varia- 
tion in the bonding arrangement of Ge and Se sites above 
x « 0.10 deviates from a chemically ordered continous 
random network (CRN). The non-random character of 
the IP implies that some form of ordering or structural 
correlation is present. In this Letter, we report a new 
electronic measurement and show that our ab initio mod- 
els of the GeSe glasses reproduce this and other features 
of the IP, and we describe the spatial correlations re- 
sponsible for the electronic signature of the IP. Our work 
shows that MD simulations can provide novel insight into 
the IP. 

In view of the experimental evidence that special spa- 
tial correlations exist in the IP, it is natural that re- 
searchers have sought out structural manifestations of the 
IP. Studies on the first sharp diffraction peak (FSDP) of 
the Ge-Se glasses indicated that the inverse peak position 
and peak area of the FSDP exhibit the telltale "flatten- 
ing" in the IP composition range[7j [8|]. However, recent 
work by Shatnawi and co-workers could not confirm these 



observations [9]. The same authors also conducted high- 
energy x-ray diffraction and extended x-ray absorption 
fine structure analyses on Ge^Sei-^ , but could not find 
structural fingerprints of the IP in the first atomic shell 

0. 

Thorpe et al. [l0(] advanced the first theoretical picture 
of the intermediate phase, in which they invoked the no- 
tion of "self-organization" - the tendency of the networks 
to develop the spatial correlations necessary to create the 
observables. Building upon earlier work on rigidity tran- 
sitions in glasses by Phillips an d Thorpe Q, they 
showed that the inclusion of extra bonding constraints 
in a floppy network gives rise to a rigid but stress-free 
bonding network, before the network transforms into a 
stressed rigid phase. Later Micoulaut and Phillips [l2| . 
by constructing networks of Ge^Sei-^ using the size in- 
creasing cluster approximation (SIC A), showed that such 
a stress free network can be thermodynamically stable. 
Chubynsky and coworkers studied rigidity percolation in 
model systems and demonstrated that a form of IP was 
possible, and potentially generic near the rigidity perco- 
lation threshold[13]. These theories have advanced our 
understanding of the IP. However, none of them offer 
a quantitative comparison with experiments, and none 
of them are based on models formed with realistic in- 
teratomic interactions. Recently, we have presented first 
principles MD models of these glasses over a wide compo- 
sition range including the IP window [16[ . In these mod- 
els, the evolution of structural parameters such as mean 
coordination and concentrations of corner-sharing (CS) 
and edge-sharing (ES) tetrahedra show a non-monotonic 
behavior. These parameters deviate from a chemically 
ordered CRN, and saturate in the composition range 
which coincides with the IP window, suggesting a link 
between experiments and the models. To a significant 
degree our work supports the picture of Micoulaut and 
Phillips [IH, inferences of Boolchand Q, and suggests the 
nature of the self-organization described in Ref. |l0|- 

In this study, we examined several compo- 
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FIG. 1: (Color online) Black points are the K edge WL po- 
sitions for Se and Ge atoms (error ±0.01eV) and red points 
are the average conduction edge energy (E c ) obtained from 
the models. 



sitions of bulk Ge^Sei-^ glasses, with x = 
0.16, 0.18, 0.19, 0.20, 0.21, 0.23, 0.24, 0.27, 0.30, 0.33. The 
glasses were prepared with a melt- quenching method 
XANES experiments were performed on beamline 9-BM- 
B of the Advanced Photon Source at Argonne National 
Laboratory. Ge K-edge [E = 11.104/ceV) and Se K-edge 
(E = 12M8keV) XANES spectra were collected in 
transmission mode. The bulk glasses were ground into 
fine powders, which were then spread onto kapton tape. 
These powders were freshly prepared just before each 
of the XANES experiment to minimize contact with 
air. Twelve layers of such tape were stacked to optimize 
the signal-to- noise ratio of the spectra. A Geo.33Seo.67 
powder sample was used as a reference to calibrate the 
absorption edge positions for all the samples. 

Strong x-ray absorption peaks are observed at Ge and 
Se K edges. These two peaks, usually called white lines 
(WLs), originate from the electronic transition between 
the core Is orbital and the lowest unoccupied electronic 
states (viz., the conduction band). Because of the fully 
screened core holes, the WL is an indication of the con- 
duction band structure [HI . Any change in the WL peak 
position implies a relative shift of the electronic states in 
the conduction band. Fig. [1] shows the variation in the 
relative K-edge WL positions for Se and Ge atoms with 
Ge content. In case of Se, the WL position decreases to a 
'plateau' for x < 0.33 with a sudden increase at x = 0.24 
and x = 0.19 close to the IP window. A similar profile 
is observed for the WL positions for Ge atoms, though 
the shift is small compared to the shift in WL position 
for Se atoms. Note that the Ge K-edge and Se K-edge 
XANES spectra were obtained independently, so the ob- 



served shift in WL positions at x = 0.19 and x = 0.24 
is not an artifact of the experiment. Since we have used 
the same reference for all the samples, the uncertainty 
in the WL positions was determined by combining the 
errors from repeated measurements of the same sample 
as well as from measurements of different samples with 
the same composition. 

Atomistic models of the IP were obtained using the ap- 
proximate ab initio density functional code FIREBALL 
[171 ]. We generated a sequence of 500-atom models with 
x = 0.15,0.18,0.22,0.23,0.25. Atoms were randomly 
placed in a cubic cell with fixed volume and experimen- 
tal density, and equilibrated above the melting point at 
1500 K for about 3.5 ps. Then they were cooled to 
400K over 4.5 ps, using velocity rescaling. Finally they 
were steepest-descent quenched to obtain relaxed con- 
formations. The models reproduce the measured pair- 
correlation functions obtained from X-ray diffraction [16]. 
Models for x = 0.10, 0.20, 0.33 are from Refs. [H, H- 

The shift in the WL positions obtained from XANES 
can be directly linked to the shift in the conduction band 
obtained from the models. In Fig. [I] red points are ob- 
tained by calculating the average conduction edge en- 
ergy (E c ) = J cb ep{e)de 1 where p(e) is the electronic den- 
sity of states and cb specifies a quadrature range starting 
at the band edge, extending 1.5 eV into the conduction 
states [20]. The absolute values are shifted by a constant 
for comparison. The similarity between the predicted 
conduction band energies and the experiments is pleas- 
ing: not only the relative shift in (E c ) is comparable to 
the experiment but it also shows a sudden increase at 
x = 0.2 and x = 0.23, a feature which is quite clear in 
the experiment as well. A small discrepancy in the posi- 
tions of the two peaks from that of experiment could be 
due to the small size and rapid quench rates used in the 
preparation of these models. Apart from that, the over- 
all profile of the shift in the absorption edge with the Ge 
content is clearly reproduced in the atomistic models. 

To elucidate the spatial character of the localized states 
at the conduction tail, the three lowest eigenstates at the 
conduction tail for compositions around the IP window 
are shown in Fig. [2j It can easily be observed that the 
reduction of the (E c ) for x > 0.10 is correlated with the 
increase in the appearance of Ge sites (green) around the 
sites contributing to the conduction tail states. Clus- 
tering of Ge creates short strained Se chains which con- 
tribute to the localization of the tail states. The effect 
has increased for x = 0.18 which corresponds to further 
lowering of (E c ). At x = 0.20 though, we see less cluster- 
ing of green sites around the eigenvectors which correlates 
with the sudden increase in (E c ) at x = 0.20. Similarly 
for x = 0.22 and 0.25, the eigenvectors are surrounded 
by more compact Ge clusters as compared to x = 0.23 
where the clusters are relatively open. Visual inspection 
shows some indication of a correlation between the shift 
in the conduction edge energy (E c ) and the clustering of 




FIG. 2: (color online) Conduction edge electron eigenstates in Ge x Se\- x . Black and red sites are Se and Ge atoms respectively 
which contribute to the lowest three eigenstates (L, LI and L2, where L stands for "lowest unoccupied molecular orbital" state) 
- states at the conduction tail, and white and green are the neighboring Se and Ge atoms. 



Ge sites. The smaller variation in Ge WL reflects the 
more uniform Ge bonding environments relative to Se. 

To explore these correlations further, we divide the cell 
in three regions. Region A consists of CS/ES tetrahedra 
(n = 1 chains), region B consists of short chains Se n 
(n — 2, 3, 4, 5) and region C is composed of longer Se 
chains representing the amorphous Se background. Fig 
[3^, shows the concentration of atomic sites belonging to 
three regions, which contribute to the valence and con- 
duction localized states. In the IP range, tail states are 
mostly localized in B region, which is understandable as 
this region appears around the clusters of CS/ES tetra- 
hedra due to cross-linking and thus becomes strained as 
the size of the clusters and the content of this region in- 
creases to a maximum in the IP window. Fig [3]o shows 
the average lengths of n = 1 and n = 2 chains, the strain 
is apparent in terms of the 'bending' of these chains. The 
average length of n = 1 chain follows a similar profile as 
the shift in the conduction edge energy (Fig. [I]). In or- 
der to see the effect of the bending of these chains on the 
electronic spectrum, we studied a toy model consisting of 
these two chains attached to each other [22| . The decrease 
in the lengths of these chains shows a clear downward 
shift in the eigenvalues above the fermi level compared 
to the relaxed system, which suggests that the antibond- 



ing states of this system are more sensitive to the bond 
angles. The two peaks in the conduction edge energy at 
x = 0.20 and x = 0.23 thus can be understood in terms 
of the variation in the lengths of these units. The role of 
filaments in tail states has recently been discussed (il). 

A comprehensive analysis of the topological evolution 
with the Ge content will be presented elsewhere [22[. In 
compact form, addition of Ge in a-Se network nucleates 
clusters of CS/ES tetrahedra (region A), which can be 
considered as harbingers of g-GeSe2 since CS/ES tetra- 
hedra are the building blocks of this network. Due to the 
cross-linking at Ge sites, these clusters are surrounded by 
short Se n chains (region B). Increase in cluster size also 
increases the concentration of Se n chains at the cost of 
the background a-Se (region C). Region B acts as a 'bar- 
rier' against the formation of GeSe2 (volumes of type A), 
thus resisting the percolation of region A. The IP appears 
when most of a-Se background "C" is transformed into 
short chains "B" and the system has many subcritical 
fragments of g-GeSe2 phases and the remainder is short 
Se n chains. Appearance of such a mixed phase is the 
cause of the apparent delay in the evolution of the net- 
work [16] [22], characterized in various experiments as a 
flattening of observables for x in the IP window. Further 
increase in the Ge content causes cross-linking in region 
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FIG. 3: a) (Color online) Concentration of sites represent- 
ing the three regions, which contribute to ten localized states 
at the valence and conduction tails for different compositions 
(see text), b) Variation of the average lengths of n = 1 and 
n — 2 chains with the Ge content. Vertical dotted lines indi- 
cate the width of the IP. 

B, and thus region A percolates all over the system. It 
seems that the evolution of the region B defines the width 
and the position of the IP on the composition scale. As 
Ge atoms tend to form clusters rather then getting dis- 
persed, it ensures the persistence of region B in the IP 
window. Clustering of Ge sites can be inferred from the 
increase in the size of FSDP with Ge content 9] which is 
mainly contributed by Ge-Ge correlations [23]. The ap- 
pearance of short Se chains is apparent from the study of 
the evolution of the Raman Se chain mode which splits 
into two modes above x = 0.19 [24]]. Higher frequency 
modes are attributed to short Se chains. In this picture 
the existence of the IP is due to the appearance of short 
Se chains. Severing these chains would result in the col- 
lapse of the IP. This may explain the reduction in the 
width of the IP due to the addition of I for the x = 0.25 
composition [25[. Also since in the IP range the band 
tail states reside mainly on short Se chains (fig. [3^i), the 
IP would also be sensitive to the light illuminations as 
micro- Raman studies suggests [26| . 

In summary, we have presented an experimental elec- 
tronic signature of the IP in Ge^Sei-^ glasses from 
XANES measurements. First principles models repro- 
duce the composition-dependent shift in the conduction 
edge. The shift is explained in terms of the change in the 
lengths of short Se n chain units. An atomic picture of 
the IP is presented with clear links to the experiments. 
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